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Physical Observables from Lattice QCD at Fixed Topology* 
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Because present Monte Carol algorithms for lattice QCD may become trapped in a given topological charge 
sector when one approaches the continuum limit, it is important to understand the effect of calculating at fixed 
topology. In this work, we show that although the restriction to a fixed topological sector becomes irrelevant in 
the infinite volume limit, it gives rise to characteristic finite size effects due to contributions from all 8- vacua. We 
calculate these effects and show how to extract physical results from numerical data obtained at fixed topology. 
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1. TOPOLOGICAL CHARGE SECTORS 
AND 9- VACUA 

In principle, one should calculate observables 
in QCD with finite quark masses in a fixed 9- 
vacuum. However in practical lattice calcula- 
tions, algorithms that change the gauge field con- 
figuration in small steps tend to become trapped 
in a fixed topological charge sector because they 
cannot overcome the action barriers between sec- 
tors. One sees concrete evidence of this prob- 
lem, for example, from the large equilibration 
times for the topological charge required in hy- 
brid Monte Carlo calculations, which increase as 
the quark mass is decreased [EJ. In addition, im- 
proved pure gauge actions, such as DBW2, re- 
duce the low eigenmodes that are undesirable for 
overlap and domain wall fermions by suppress- 
ing small instantons and thereby reduce tunneling 
between sectors. Since we know that the restric- 
tion to fixed topology becomes irrelevant in the 
infinite volume limit, we have calculated the the 
finite volume dependence of this restriction and 
thereby show how to extract physical results from 
practical calculations at fixed topology. 
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The partition function of the QCD Hamilto- 
nian H with eigenstates \n, 9) in a given sector 
may be written in terns of the lowest eigenvalue 
Eq(9) = Veo(9) in the large volume, low temper- 
ature limit as follows: 
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The partition function at fixed topological charge 
Q is an integral over all ^-vacuum sectors. Ex- 
panding the ground state energy 



eo(9) = e (0) 
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and performing a saddle-point expansion to 
0(l/p 2 V 2 ), 
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Using Z e = JVAVVVy e -siA,^]-ieQ[A]^ the 
parameters in the expansion of eo(9) are 
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The topological susceptibility, \t, an d fourth 
moment, 7, are parameters characterizing the 
vacuum that can be calculated in lattice QCD. 
Once they are known, we may use eq. (|l|) to cal- 
culate any observable at fixed Q and obtain finite 
volume corrections. 

2. MASS SPECTRUM IN A FIXED 
TOPOLOGICAL SECTOR 

As a concrete example, consider the calculation 
of hadron masses from the large time behavior 
of two-point correlation functions. At sufficiently 
low temperature, the two-point correlation func- 
tion at fixed Q of operators O with the appropri- 
ate quantum numbers to create the hadron states 
of interest may be written 

(0(h)0(t 2 )} Q ' 



d Q e i9Q~f3Ve {0) 



x \{Q,e\0\l,Q)\ 2 e- M ^-^\ 

where M(9) is the lowest hadron mass. Perform- 
ing a saddle-point expansion, the effective mass 
used to measure the hadron spectrum is 

-^log(<0(ti)0(ta)>g) 
fd0e-^Q 2 ^ 2+ie Qf{8) M{6) 
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where f(9) = \(0,9\O\l,9)\ 2 e~ M ^^- t2 ^ 
contains all the time dependence. Expanding 
M{9) = M(0) + \M"{Q)9 2 + ^M""(0)6> 4 + • • • 
and f(6) = /(0) + ±/"(0)# 2 + i/""(0)^ 4 4 
obtain the desired expansion, 
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The moments required in eq. (|3|) are 
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where (Q 2 ) = /3Ux* by eq. (§). Note that to 
leading order, all dependence of the mass shift on 
f(0) exactly cancels and therefore the result is 
independent of t\ — t% ■ 

Our leading order result for the Q-dependence 
of the mass is thus 

1 M"(0) I" Q 2 

^~Wx~t [ ~Wxt 

Therefore, if a calculation in a volume j3V is 
trapped in a Q sector, the error is of order 
1/((3V). Furthermore, one can measure Mq in 
several sectors and at several space-time volumes 
and thereby extract Af(0) , M"(0), and \t us- 
ing eq. (0) . Also note that when averaged over Q 
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with the distribution e 2 <- Q2 '> . 
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3. 9 DEPENDENCE IN THE CHIRAL 
AND LARGE N c LIMITS 

It is possible to make some general argu- 
ments concerning the ^-dependence in QCD. The 
transformation ip' — exp( i 2 ~wjl5)' l l } shifts the 9- 
dependence to the mass term, 

to** — > * {mcos(9/Nf) + imsin(9/Nf)"f 5 } * 

and the contribution to the fermionic measure 
cancels the FF-term. By parity, the sin(9/Nf) 
term only contributes in even orders, so at lin- 
ear order in the quark mass, the only effect of 9 
enters through the change of the mass term to 
mcos(9/Nf). Hence, for example, to leading or- 
der in the chiral limit, the pion, nucleon and 7/ 
masses must have the forms: 

M 2 (9) = M 2 (0) + M 2 (0)[cos(9/N f )-1] 

M N (9) = M N (0) + Cl M 2 {0)[cos(9/N f )-l} 

+ c 2 M3(O)[cos 3 / 2 (0/A / )-l] 

M 2 {9) = M 2 , (0) + hM 2 (0)[cos(9/N f ) - 1] 

from which the Q-dependence follows directly 
from eq. (f|). In the large iV c limit, it can be 
shown that b\ = 1 so that the n and rf have 
the same ^-dependence and thus the same Q- 
dependence. 

Using chiral perturbation theory, one can ex- 
plore the chiral and large N c limits more system- 
atically. Using the effective Lagrangian proposed 
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by Witten j| to treat the Ua(1) anomaly, we ob- 
tain the following results for the tt and rj masses: 
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where 

e (6) = -m(W) cos ^ - ^) + $rf 2 . 
Solving explicitly for t/q , we also note that in the 
limit m — > 0, we recover cos (^~) an< 4 m the 

> 1 as 



limit N c — > oo we obtain cos ( 
AT C — > oo as expected. 
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4. Q DEPENDENCE IN AN 
INSTANTON GAS 

Recently, the dependence of the it and rj 
masses on Q has been studied in lattice calcu- 
lations JjJ, which for practical reasons, are suffi- 
ciently far from the chiral limit that the previous 
results from chiral perturbation theory are inap- 
plicable. Each mass was evaluated in two sectors, 
one with \Q\ < 1.5 and the other with \Q\ > 1.5. 
Whereas the pion mass in both sectors agreed 
within statistics of a few percent, the rj mass was 
of the order of 15 % heavier in the \Q\ > 1.5 sec- 
tor than in the \Q\ < 1.5 sector. Motivated by the 
success of the Veneziano- Witten formula || [Q re- 
lating the rj mass to fluctuations in the topologi- 
cal charge, we obtain a qualitative understanding 
of the Q dependence from the fluctuations in the 
topological charge arising in an instanton gas. 

The vertex generating the shift in the rj mass is 
proportional to the number of instantons, N, plus 
the number of antiinstantons, N. Assuming inde- 
pendent Poisson distributions with (JV) = (N) = 
A, the probability of having N instantons and N 
antiinstantons with the constraint N — N — Q is: 
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Summing over N and N and distinguishing A and 
A' for subsequent convenience, we write the gen- 
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Zq(X,X') 
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Differentiation of lnZ<g(A, A') with respect to A', 
setting A = A', and noting 2A = (Q 2 ) — XtftV 
yields the desired result for the density of instan- 
tons plus antiinstantons at fixed Q: 
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Eq. ^| is a simple and physically appealing re- 
sult. Since M 2 = M% + p and /i oc {N +y )Q , 
which equals Xt when averaged over the distri- 
bution P(Q) 
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e , the rj mass is con- 
sistent with the Veneziano- Witten formula. The 
finite volume corrections provide the desired Q 
dependence. 

To compare with the lattice results of ref |Q, 
we note that (3V = 6.81/m 4 , X t = 0.70/m~ 4 , 
and (Q 2 ) = PV\t = 4.75. In the chiral limit, 



M*ff(Q) = M^O) 

the shift at Q 2 = is - NU / (AN 2 (Q 2 )) = 
0.013^/^. Hence, the effect of calculation at fixed 
topology is of the order of 1%, consistent with the 
lattice 7r results. The analogous chiral estimate 
for the rj is of order (M^/M^i) 2 x 1%, in strong 
disagreement with lattice results. However, us- 
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is 



ing the instanton gas result, / dqP(Q) 

0. 65 for Q < 1.5 and 0.75 for Q > 1.5, yielding 
5M rj i /M n > — 8%, which is of the order of magni- 
tude of the observed Q dependence. 
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